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We search for a new light gauge boson, a dark photon, with the DO experiment. In the model we 
consider, supersymmetric partners are pair produced and cascade to lightest neutralinos that can 
decay into the hidden sector state plus either a photon or a dark photon. The dark photon decays 
through its mixing with a photon into fermion pairs. We therefore investigate a previously unex- 
plored final state that contains a photon, two spatially close leptons, and large missing transverse 
energy. We do not observe any evidence for dark photons and set a limit on their production. 

PACS numbers: 95.35. +d, 12.60.Jv, 14.80.Ly 



Hidden valley models [T] introduce a new hidden sec- An important subset of hidden valley models also con- 
tor, which is very weakly coupled to the standard model tain supersymmetry (SUSY), a fundamental symmetry 
(SM) particles, and therefore can easily escape detection. between fermions and bosons postulating the existence 
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of SUSY partners. At colliders, in the case of i?-parity 
conservation [2] , superpartners are produced in pairs and 
decay to the SM particles and the lightest superpartner 
(LSP). The LSP is a stable, weakly interacting particle, 
and can not be detected in collider detectors. 

Recently, these models were called upon to explain the 
results of several cosmic ray detection experiments |4] . 
Taken together with other experiments, including new 
results from Fermi/LAT [5], there is evidence of an ex- 
cess of high energy positrons and no excessive produc- 
tion of anti-protons or photons. The excess can be at- 
tributed [6] to the dark matter particles annihilating into 
pairs of new light gauge bosons, dark photons, which 
are force carriers in the hidden sector. The dark pho- 
ton mass can not be much larger than 1 GeV to give 
rise to Sommerfeld enhancement [7] of the dark matter 
annihilation cross section, and not to decay into neu- 
tral pions and/or baryons. The masses of the hidden 
sector states are also around 1 GeV, with mass split- 
ting around MeV, thus providing a possible explana- 
tion of the DAM A [5] signal through "inelastic Dark 
Matter" scenarios. Dark photons decay through mix- 
ing with photons into SM fermions with branching frac- 
tions that can be calculated from the measurements [3] 
of R = a(e + e~ — > hadrons)/er(e + e~ — ^ fi + fi~), and 
strongly depend on the dark photon mass. For dark 
photon masses (m 7D ) below the dimuon threshold of 
~ 200 MeV, only decays into electrons are possible. For 
m lD ~ 0.5 GeV the decay rates into electrons and muons 
are approximately 40% each. The lowest value of the lep- 
tonic branching (3.7%) occurs if the dark photon mass is 
accidentally equal to that of the <fr meson. 

In this Letter we will follow the phenomenological sce- 
nario developed in [10 . A diagram of a possible process 
at the Fermilab Tevatron Collider is shown in Figure [T] 
Gauginos are pair produced and decay into SM parti- 
cles and the lightest neutral gaugino (neutralino, Xi), 
which in turn decays with comparable branching ratios 
into either a hidden sector darkino X (which is the LSP), 
and a photon, or into darkino and a dark photon (713). 
Hadronic dark photon decays are overwhelmed by SM jet 
backgrounds. Thus, we only consider dark photon de- 
cays into isolated electron or muon pairs. Both darkinos 
escape detection and result in large missing transverse 
energy The branching fraction of the neutralino 

into the dark photon, B — Br(xi —> IdX), is a free pa- 
rameter of the model. If it is small, the decays into a 
photon dominate, and signature is the same as of GMSB 
SUSY jTT] with the neutralino as next-to-lightest super- 
partner (NLSP). Larger values of B give rise to events 
where one of the two neutralinos decays into a dark pho- 
ton, resulting in a final state with one photon, two spa- 
tially close (and therefore not satisfying traditional iso- 
lation requirements) leptons and large $ T . 

This Letter describes a search for this, so far unex- 
plored, final state in pp collisions at a center of mass 




FIG. 1: One of the diagrams giving rise to the events with 
a photon, dark photon (70), and large missing energy due to 
escaping darkinos (X) sit the Fermilab Tevatron Collider. 

energy of 1.96 TeV recorded by the DO detector [12] 
at the Fermilab Tevatron Collider. As is described be- 
low, our search is optimized for low dark photon masses, 
m lD < 2.5 GeV. Note that it is also sensitive to the case 
where the neutralino decays into a hidden state Y with 
somewhat higher mass than the dark photon. The Y may 
cascade down to the darkino through other hidden states 
which may be long-lived and can result in the emission of 
highly collimated low energy SM particles, some of which 
could be leptons. Most of the energy of the Y will stay in 
the hidden sector and therefore the high $ T should not 
be substantially reduced. This analysis is also sensitive 
to another possible scenario, proposed in |13j . in which a 
light axion that decays into muon pairs takes the place 
of the dark photon in the decays described above. 

Data for this analysis correspond to an integrated lu- 
minosity of 4.1 fb _1 after application of data quality 
and trigger requirements. Events must satisfy a set of 
high transverse energy (Ex), single electromagnetic (EM) 
cluster triggers which are fully efficient for photons with 
E T > 30 GeV. 

EM clusters are selected from calorimeter clusters, 
built using the simple cone algorithm of radius 1Z = 
V(Ar?) 2 + (A0) 2 = 0.4 121, by requiring that the frac- 
tion of the energy deposited in the EM section of the 
calorimeter, EMy rac , is above 95% and the calorimeter 
isolation variable I = [E tot (0A) - E EM (0.2)]/E EM (0.2) 
is less than 0.2, where E to t(0A) is the total energy in 
a cone of radius 1Z — 0.4, corrected for the underlying 
event contribution, and E em (0-2) is the EM energy in 
a cone of radius 1Z — 0.2, which is taken to be the EM 
cluster energy. 

Photon candidates are selected from central calorime- 
ter (| 77 1 < 1.1) EM clusters which have (i) EM /rac > 97%, 
(ii) X < 0.07, (iii) a shower shape consistent with that 
of a photon, and (iv) the scalar sum of the transverse 
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momenta (pr) of all tracks originating from the primary 
vertex in an annulus 0.05 < 1Z < 0.4 around the clus- 
ter less than 2 GeV. Additionally, we require that pho- 
ton candidates are not spatially matched to activity in 
the tracker. The tracker activity can either be a recon- 
structed charged particle's track or a density of hits in 
the silicon microstrip and central fiber trackers consistent 
with a track. The EM clusters that do not have matched 
activity in the tracker, but fail other photon selection 
criteria, are dominated by jets that have fragmented into 
neutral pions, and are referred to below as fake photons. 

We search for dark photons in events with at least one 
photon with E T > 30 GeV and $ T > 20 GeV {$ T 
is computed using all calorimeter cells and corrected for 
EM and jet energy scales). Dark photon candidates are 
formed by selecting pairs of oppositely charged spatially 
close (TZ < 0.2) tracks that originate from the same 
point (|Az| < 2 cm) along the beam line. The lead- 
ing (trailing) track p T is required to exceed 10 (5) GeV. 
We then require the scalar sum of pr of all tracks ex- 
cluding the pair in a cone of radius 0.4 centered on the 
pair momentum direction to be less than 2 GeV. To re- 
duce the QCD background we require that each track 
must have its azimuthal angle not aligned with a photon, 
0.4 < A(f> J trac ), < 2.74. In rare cases, when there is more 
than one such pair in the event, we select the one with 
the highest trailing track px- 

For a dark photon decaying into a pair of electrons, 
the calorimeter depositions overlap, so we require that 
the dark photon candidate matches an EM cluster with 
E T > 10 GeV, EM /rac > 97%, and 1 < 0.1. For 
the dimuon decay mode, we require that at least one 
of the tracks is matched to a reconstructed muon, and 
the energy deposited in the calorimeter in the annulus 
0.1 < TZ < 0.4 is below 3 GeV. 

Dark photons would manifest themselves as a narrow 
peak in the lepton pair invariant mass distribution. We 
use a Monte Carlo simulation to characterize the mass 
resolution, as well as the efficiency to reconstruct the 
events. SUSYHIT [T3] is used to calculate masses and de- 
cay probabilities for the GMSB SUSY [IB] model, known 
as Snowmass Slope SPS 8, and produce the Les Houches 
Accord [T7] card files. These files are modified to intro- 
duce neutralino decays to a dark photon. Events with 
one of the two neutralinos decaying into a dark photon 
and the other decaying into a photon are generated with 
PYTHIA [18] using CTEQ6L1 parton distributions [19] and 
are passed through the full GEANT-based [20] detector 
simulation and the same reconstruction chain as the data. 
Following [21 , the leading order (LO) signal cross sec- 
tions calculated by PYTHIA are scaled to match the next- 
to- leading order (NLO) prediction using fc-factor values 
extracted from [22] . The event kinematics depends on 
the mass of the dark photon and the masses of super- 
partners, resulting in slight variations in mass resolution 
and selection efficiency. Typical values are 5% and 12%, 



respectively. 

There are three types of SM processes that contribute 
to our data sample: 

Bl QCD events with real or fake photons and mismea- 
sured $ T . These contain jets or photon conversions 
faking the dark photon. 

B2 W — > e/fi v plus a real or fake photon. The dark 
photon is faked by a accidental overlap of a high 
Pt track with the lepton. 

B3 W — ► tv — > 3h v plus real or fake photon. One 
of the particles from t lepton decays is lost or very 
soft, and the remaining t decay products fake the 
dark photon. 

We study dark photon candidate mass distributions in 
three control samples where we do not expect dark pho- 
tons to appear. The QCD^f control sample is selected 
by reversing the $ T cut. The QCD_jet sample is se- 
lected by using the same criteria as the QCD_j sample, 
but requiring a fake photon instead of a photon. Fi- 
nally, the QCD_W sample requires a fake photon and 
$ T > 20 GeV. All three have contributions from Bl, 
although the relative fraction of multijets, single photon 
production, and diphoton production varies among the 
three control samples. Backgrounds B2 and B3, how- 
ever, can only significantly contribute to the QCD_W 
sample. We observe no difference between the dark pho- 
ton candidate mass distributions in the three control 
samples. We therefore conclude that the background to 
dark photon production is dominated by Bl, and use the 
average shape of the dark photon candidates mass distri- 
butions in all control samples as our background model. 

The dark photon candidate invariant mass distribu- 
tions in the signal sample are shown in Figure [2] sepa- 
rately for the electron and muon channels, together with 
the expected contribution from dark photons with a mass 
of 1.4 GeV. 

We see no evidence of a dark photon signal and pro- 
ceed to set limits on its production. To set limits we use 
the standard DO likelihood fitter [23] that incorporates 
a log-likelihood ratio (LLR) statistic method [24 . The 
value of CL S is defined as CL S = CL s+ b/CLi,, where 
CL s+ b and CL\> are the confidence levels for the sig- 
nal plus background hypothesis and the background-only 
(null) hypothesis, respectively. These confidence levels 
are evaluated by integrating the corresponding LLR dis- 
tribution populated by simulating outcomes via Poisson 
statistics. Systematic uncertainties are treated as uncer- 
tainties on the expected number of signal and background 
events, not the outcomes of the limit calculations. This 
approach ensures that the uncertainties and their corre- 
lations are propagated to the outcome with their proper 
weights. The limit is set by simultaneously fitting dilep- 
ton invariant mass distributions in data for the muon and 



6 



^ 35i_ 1 1 1 1 i 1 1 1 1 i 1 
u : D04.1fl 

O 30r m+ii . 




> 

O 



2.5 



1.5 2 2.5 
(GeV) 



FIG. 2: Observed mass distributions in the signal region are 
represented as points with error bars, the background esti- 
mation is shown as filled band, and an example signal for 
m 7D =1.4 GeV plus background is shown as the solid his- 
togram for the dimuon channel (a) and the dielectron channel 
(b). 



electron channels to the signal and background predic- 
tions for each signal point, defined by the dark photon 
and the lightest chargino masses. For each dark pho- 
ton mass the background is normalized outside of the 
expected signal region. The systematic uncertainty on 
the signal reconstruction efficiency (25%) is dominated 
by the uncertainty to reconstruct the two spatially close 
tracks from the dark photon decays (20%). The latter 
was cross-checked with data using tau decays and con- 
verted photons in radiative Z — > fifjfy decays. We also 
took into account the uncertainty on the total integrated 
luminosity (6.1%) and the effect of varying the dark pho- 
ton mass resolution by 10%. 

We interpret the cross section limits as limits on the 
lightest chargino mass as a function of the dark photon 
mass and the neutralino branching fraction. For B = 0.5 
the excluded region of chargino and dark photon masses 
is shown in Figure [3] In Figure [4] we display the chargino 
mass limit as a function of B for three representative 
dark photon masses: 0.2 GeV (only the electron channel 
is open), 0.782 GeV (low branching fraction into leptons 
due to u) and p mesons), and 1.5 GeV. Our previous limit 
on the GMSB SUSY in the diphoton final state [21] is 
directly applicable to the model considered in this Let- 
ter, although it does not probe the dark photon mass. 
The corresponding exclusion contours are shown in Fig- 
ures |3] and 3J 

To summarize, we search for a previously unexplored 
final state consisting of a photon, two spatially close lep- 
tons from hypothetical dark photon decays and large 
missing energy. We find no evidence for such events, and 
set limits on their production in a benchmark model [TO] ■ 
For dark photon masses of 0.2, 0.782, and 1.5 GeV we 
exclude chargino masses of 230, 142, and 200 GeV, re- 
spectively. 

We would like to thank Scott Thomas and David Shih 
for many inspiring discussions and help with the signal 
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FIG. 3: The excluded region of possible masses of the light- 
est chargino and the dark photon for B = 0.5 are shown as 
the shaded region. The expected limit is illustrated as the 
dash-dotted line. The vertical black line corresponds to the 
exclusion from the diphoton search ,21.. 
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FIG. 4: The dependence of the excluded chargino masses 
on the branching ratio of the neutralino into a photon are 
given for dark photon masses of 0.2, 0.782, and 1.5 GeV. The 
black contour corresponds to the exclusion from the diphoton 
search |21| . 
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